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Abstract 



m 

J> ' We report on the result for the charm quark mass as obtained from our lat- 

tice QCD computation in the quenched approximation. Our result is m^r (m c ) = 
1.26(4)(12) GeV. 



PACS: 11. 15. Ha (Lattice gauge theory), 12.38.Gc (Lattice QCD calculations), 14.65. Dw (Charmed quark) 



a,: 

Ph! 1 Introduction 

1.1 Status of the charm quark mass calculations 

/^ ■ During the last few years a great effort has been devoted to the precise determination of 

the quark masses. Recent reviews about the present situation concerning the computation 
of the light quark masses can be found in refs. |l|, @ . As for the heavies, most of the studies 
done so far were focused on determining the 6-quark mass. The charmed quark escapes the 
precision computation mainly because it is too heavy for the chiral perturbation theory to 
apply, and yet too light for the heavy quark expansion to set in. 

A complete account of the presently available estimates of the charm quark mass value 
is given in the PDG review ||. They estimate the charm quark mass to be in the range [] 

1.15 GeV < mfS( mc ) < 1.35 GeV . (1) 

Very recently, two new QCD sum rule computations of this quantity appeared 0, ^| . After 
improving the calculation of the moment QCD sum rules for the charmonium states, the 
new result of ref. [[J is m^lmc) = 1.23(9) GeV. Adopting quite a different QCD sum rule 
technique, in ref. 0] m^ s (m c ) = 1.37(9) GeV was obtained. 

On the side of the lattice QCD computations there were a few attempts to compute the 
charm quark. 

1 Ref. [Bj also contains a complete list of references. 



By combining the QCD sum rule methodology with the lattice QCD computation of 
the moments of the heavy-heavy vector current correlation function, the authors of 
ref. || obtained m}f s (m c ) = 1.22(5) GeV. More complete discussion of the systematic 
uncertainties involved in their computation has been made in ref. 0. Unfortunately, 
this (elegant) method has not been followed by the other lattice groups. 

From the lattice QCD with Wilson quarks and by using the Ward identities, in ref. || 
a much larger value has been obtained, namely m^(mc) = 1.71(3)(20) GeV. If we 
use their recent re-estimate of the mass renormalization constant 0, the above value 
gets down to mf(m c ) = 1.50(3) (18) GeV. 



Finally, in ref. [[TIJ the charm quark mass has been computed by adopting the Fermilab 



way of interpreting the lattice QCD beyond the lattice cut-off. The quoted result is 
m ^S(m c ) = 1.33(8) GeV. 

Since the actual situation with the lattice results is not clear, we decided to take advantage 
of the presently available non-perturbatively determined parameters, which are necessary 
for the complete elimination of the discretization effects that are linear in the lattice spacing 
(i.e. of 0(a)), and to make an estimate of the charm quark mass from the data that we 
obtained by working with the (relatively) fine grained lattice (a ~ 0.07 fm). Our result is 

m fS(m c ) = 1.26(3) (12) GeV . (2) 

1.2 Computation of the charm quark mass on the lattice 

To compute the charm quark mass, we rely on the standard two strategies: 
• Vector Ward identity (VWI) 



m Q 



(ji) = Z m ({i)m ( Q W1) = Zjjjk)- ( I , (3) 

6 V K Q K cr J 



where the critical value of the Wilson hopping parameter (« cr ) corresponds to the 
chiral limit. Q Since the lattice cutoff is finite and the charm quark mass is not 
negligible, it is highly important to improve the renormalization constant out of the 
chiral limit: 

Z m (^ = -^—(l + b m m^ W1) ), (4) 

z s (N 



where the value of b m = — 6^/2, has been determined non-perturbatively in ref. [11 
and Z s (/i) is the renormalization constant of the scalar density computed in the 
chiral limit, which will be given in the next section. 



2 Unless the physical units are explicitely displayed, all the quark masses mentioned in this section are 
assumed to be in lattice units. 



Axial-vector Ward identity (AWI) 

(J>4(z)^(o)> 
<*M h- < WI) (.) = 2-M ( ^p (x)P t (0)) • < 5 > 

where A^ = 0tm7s9 an d -P — Q75<? are the axial vector current and the pseudoscalar 
density, respectively. For the full elimination of 0(a) effects, the bare lattice axial 
current is improved in the chiral limit as A 1 = A^+caO^P, with ca already determined 
non-perturbatively fi~2| , [I~3|j . We used the symmetric definition of the derivative, e.g. 
d^f = (/(t + 1) — fit — l))/2. The mass renormalization constant is equal to Z m (fi) = 
Za/Zp(ii) and it is improved as 



r(0) / m (VWI) (VWI). 

ZrM = -Sfc- ( 1 + (6x - 6p) q t 9 ) • (6) 



Again, the value of (pA — bp) has been computed non-perturbatively [D], |T2|, whereas 
the value of Z m (//) = Z^/Zp (//), will be provided in the next section. 

2 Lattice details and results 

2.1 A few details on the lattice computation 

We work with two sets of lattice data, generated at (3 = 6.2 by using the non-perturbatively 
improved Wilson action (csw — 1-614 [H|). Each set contains 200 independent 577(3) gauge 



field configurations. The value of the critical parameter, K cr , is fixed by requiring the bare 
m q — > 0. The values of n cr , along with the values of the bare quark masses for the 
degenerate light quark combinations, as well as for the non-degenerate (heavy-light) ones, 
are given in tab. |l[ Each m*- AWI ' and m' VWI ) are obtained as 

£>4(s)pt( )) 



(7) 



™(AWI) _ 


1 X 




2 £>(z)P'(0)) 


m (VWI) = 


= i(-+---) 

4 \Kq K q K cr J 



where in the improved axial current, A\ = Q7475<? + CAd^Qj^q, we used ca = —0.04(1) |T2]. 
p~3|] . The statistical errors quoted in this work are obtained by using the standard jackknife 
procedure (with 5 configurations per jack). The quark masses listed in tab. [I] are the bare 
lattice ones, that we now need to renormalize. We next discuss the computation of the 
mass renormalization constants. 



Lattice (f3 = 6.2) 


24 3 x 64 (n cr = 0.13583(5)) 


24 3 x 48 (n cr = 0.13577(3)) 


«l-«2 


m (AWI) 


m (VWI) 


M P 


m (AWI) 


m (VWI) 


M P 


0.1344-0.1344 
0.1349-0.1349 
0.1352-0.1352 


0.0393(15) 

0.0249(10) 

0.0162(9) 


0.0391(13) 
0.0253(13) 
0.0171(13) 


0.305(2) 
0.244(2) 
0.200(3) 


0.0399(14) 
0.0241(12) 
0.0155(11) 


0.0376(9) 
0.0237(9) 
0.0154(9) 


0.307(2) 
0.245(2) 
0.200(2) 


0.125-0.1344 
0.125-0.1349 
0.125-0.1352 


0.176(2) 
0.167(2) 
0.163(3) 


0.179(1) 
0.172(1) 
0.168(1) 


0.690(1) 
0.672(2) 
0.661(2) 


0.172(2) 
0.162(2) 
0.158(3) 


0.177(1) 
0.170(1) 
0.166(1) 


0.693(2) 
0.675(2) 
0.663(2) 


0.122-0.1344 
0.122-0.1349 
0.122-0.1352 


0.221(3) 
0.212(3) 
0.207(3) 


0.228(1) 
0.221(1) 
0.217(1) 


0.786(1) 
0.768(2) 
0.757(2) 


0.217(3) 
0.209(3) 
0.202(3) 


0.226(1) 
0.220(1) 
0.215(1) 


0.789(2) 
0.771(2) 
0.761(3) 


0.119-0.1344 
0.119-0.1349 
0.119-0.1352 


0.267(3) 
0.258(3) 
0.251(3) 


0.280(1) 
0.273(1) 
0.269(1) 


0.876(1) 
0.858(2) 
0.847(3) 


0.261(3) 
0.252(3) 
0.248(4) 


0.278(1) 
0.271(1) 
0.267(1) 


0.878(2) 
0.861(3) 
0.851(4) 



Table 1: Bare quark masses obtained by using the vector and the axial Ward identity on the 
lattice (see eq. (ffi)). The results are given in lattice units. On the subset of the configurations 
gathered on the lattice 24 3 x 64, in ref. jl^/ we computed the strange quark mass. 



2.2 Mass renormalization constants 

To evaluate the mass renormalization constants in the chiral limit we use the non-perturba- 
tive method (JTJJ which allows one to compute these constants in the continuum RI/MOM 
scheme. The use of this method out of the chiral limit, however, requires the improvement 
of the off-shell quantities for which several new counterterms appear, each with a coefficient 
that is to be fixed non-perturbatively (for more details, see ref. |{L6|| ). For that reason, in 
this paper, we will employ the method to compute the renormalization constants in the 
chiral limit only [IT]. To improve the renormalization constants out of the chiral limit, we 



will rely on the results of ref. [11], [12[ in which the corresponding coefficients were computed 



non-perturbatively with the results, b m = —0.69(1), b& — b P = 0.04, at j3 = 6.2. 

We first focus on the mass renormalization constant, Zm (fJ>) = 1/Z S \jx). To that 
end, we combine the quark propagators of 4 different values of n q in 10 different combina- 
tions (degenerate and non-degenerate in light quark mass), compute the amputated vertex 
function Fs(ni, Kj', cl/j,) and impose the standard RI/MOM renormalization condition (] 



lim * r = lim rsfaw) 



Z { s\fj) Ki,a-KerZ5(«i,K 2 ;ju) Ki,a-«er zl /2 (Ki;p) Z¥ 2 (Kj;p) 



(8) 



where Z q (fi) is the quark field renormalization constant which is easily extracted by impos- 
ing the vector Ward identity on the quark propagator. Such an obtained Z m (fi) = 
1/Z S (/i) is then converted to the renormalization group invariant (RGI) form, Z^ 1 = 
Z m (/x)/c RI ' MOM (/i), by using the available mass anomalous dimension coefficients up 
to 4- loops encoded in the function c RI/,M0M (/i), computed in the same RI/MOM scheme [IS 



For consistency, we also use the 4-loop expression for the running coupling [nj|], and set 
n F = with (quenched) Aqcd = 0.25 GeV [F2(], F2TJ. A typical situation is shown in 



fig. |l](left). For every mass combination, we then extract Z^(fi) by fitting to a constant on 
the plateau 1.0 < (afi) 2 < 1.8. After extrapolating in am q = (m\ + m 2 )/2 to the 
chiral limit, we get 

Zj® 1 * = 3.391(22) . (9) 

For an easier comparison with the other determinations of this renormalization constant 
we express this result also in the MS scheme: 

Z^ m {2 GeV) = 1.332(9) . (10) 

As it can be observed from fig. |](right), the extrapolation to the chiral limit in this case is 
very smooth. 

Next, we discuss the computation of the second mass renormalization constant, i.e. 
the one needed to compute the quark mass by using the axial Ward identity (Z m (/i) = 
Z^jZp (//)). We use the proposal of ref. which, by judiciously combining the Ward iden- 



tities, allows one to alleviate the problem of the contamination by the Goldstone boson [22 



3 Besides the values of the Wilson hopping parameters corresponding to the light quark masses that we 
already listed in tab ffl, for this computation we also use the quark propagator with n q = 0.1333. 
4 The explicit form of the function c RI / MOM (/Lt) is also given in Appendix of the present letter. 
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Figure 1: Mass renormalization constant: In the left figure we show the renormalization constants 
divided by the perturbative 4-loop anomalous dimension for the case of VWI (Z m = l/Zg) and 
AWI (Z m = Zj^jZp) for a specific combination of K\ = 0.1349 and k 2 = 0.1352. On the right 
figure we show the extrapolation of the mass renormalization constants to the chiral limit. Dashed 
line depicts the quadratic extrapolation to Z 
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The renormalization condition, by which this task is achieved for the mass renormalization 
constant, can be simply written as 



7 (0) 
7 (0) 



(ji) = lim — («i, k 2 ; n) = lim 

Kl,2— »Kcr Zp Kl,2— «cr 



(VWI) 

ml 



T p (ki;p) -m 2 vwl) T P (K 2 ;p) 



[m 



(VWI) 

i 



m 



(VWI) - 



Ta{ki,K2)v) 



[ID 



p*=H* 



Obviously, we can use only the non-degenerate quark mass combinations which, for our 4 
values of K q , means 6 combinations. As in the previous case, we convert our result from 
the RI/MOM to the renormalization group invariant constant (at 4-loop level) and fit in 
the same window as before, 1.0 < (yua) 2 < 1.8 (see fig. |l](left)). With such extracted 
values for Z m , for each combination of the Wilson hopping parameters («i, k 2 ), we then 
extrapolate to the chiral limit. This is also illustrated in fig. |](right). Contrary to the 

first case, the mass dependence of the Z m is more pronounced. We extrapolate to the 
chiral limit linearly (filled square in fig. [l]) to get our central value. In addition, we perform 
the quadratic extrapolation (the result of which is depicted by an empty square in fig. |l|), 
and the difference between this and the central value is incorporated in the systematic 
uncertainty. Our result is 



^(0)rgi 



3.303(26) 



+0.000 
-0.051 ' 



which in the MS scheme reads 



Z^ )MS (2 GeV) 



1.297(10) 



+0.000 
-0.020 • 



(12) 



(13) 



This result agrees well with the one of ref. [EI], Z m (2 GeV) 



1.316(14)(17). 



2.3 Putting it all together 

Now we combine all the results from tab. p] with the renormalization constants discussed 
in the previous section to get the renormalization group invariant quark masses. At this 
stage (after including the renormalization constants) one can identify the results of the axial 
Ward identity as the sum of the heavy and the light quark mass. Since we also computed 
the light quark mass separately, we can now simply subtract it from the sum and work with 
the heavies only. The results are presented in tab. 0. In the same table we give the values 
of the heavy-light mesons for which the light quark mass has been interpolated to the light 
s-quark in a usual way (see e.g. ref. |[23|| ). The mass of the D s meson in lattice units is 
Md s ~ 0.73(3), where we use a~ l (mK*) = 2.7(1) GeV. Thus the charm quark mass is to be 
found through an interpolation between the results for kq = 0.125 and k,q = 0.122. Notice 



o 


24 3 x 64 


24 3 x 48 


«Q 


mf [AWI] 


mf [VWI] 


M Ps 


mf [AWI] 


mf [VWI] 


M Ps 


0.125 
0.122 
0.119 


1.032(13) 
1.334(16) 
1.638(20) 


0.868(3) 
1.013(2) 
1.141(1) 


0.676(3) 
0.771(3) 
0.861(3) 


1.000(14) 
1.306(16) 
1.592(19) 


0.864(2) 
1.011(1) 

1.139(1) 


0.681(5) 
0.777(4) 
0.866(4) 



Table 2: Renormalized heavy quark masses directly accessed from our lattice. 

that in tab. |2| we give also the results for a larger quark mass corresponding to kq = 0.119. 
Although this value is not necessary for our final result for the charm quark mass, it will 
be helpful in assessing the amount of the systematic uncertainties which will be discussed 
in the next subsection. 

Now, to get the value of the charm quark mass, we need to interpolate in the heavy 
meson mass to Mp> a and then simply read off the charm quark mass. To do so we need to 
choose an interpolating formula. We consider the following ones: 

(i) M Ps = a + aiUiQ + a 2 m 2 Q \ 

(ii) Mp 3 =b + h/rriQ + b 2 /m 2 Q ; 

(iii) M Ps = c + CirriQ + c 2 /m Q . 

The first (i) is the naive linear interpolation (a 2 = 0), the second (ii) comes from the heavy 
quark expansion and the third (iii) is the hybrid of the two. For (iii) we obviously need 
at least three points, and thus we have to use also the heaviest of our quarks from tab. EL 
The complete situation is presented in tab. §. Since the heavier quark is more prone to 
the O ((am) 2 ) artifacts, we first concentrate our discussion on the results of the first two 
interpolations in which we set a 2 = b 2 = 0. The results of the two interpolation formulae 
are totally consistent with each other, and we choose to quote the first one as our central 
number. 



kq 



Form (i) (02 = 0) 
Form (ii)(&2 = 0) 



Form (i) 
Form (ii) 
Form (iii) 



24 3 x 64 (Lattice-I) 



m 



rgi 



[AWI] 



1.190(84) 
1.171(83) 



1.186(84) 
1.191(82) 
1.187(84) 



mfi [VWI] 



0.944(39) 
0.939(39) 



0.945(39) 
0.947(39) 
0.954(39) 



24 3 x 48 (Lattice-II) 



mjf [AWI] 



1.136(100) 

1.126(98) 



1.135(100) 

1.143(98) 

1.137(100) 



mf [VWI] 



0.929(47) 
0.923(46) 



0.930(47) 
0.932(46) 
0.931(47) 



Table 3: Charm quark mass in lattice units as obtained by using the interpolating formulae (i), 
(ii) and (Hi), as discussed in the text. The first two are obtained without using the heaviest of 
our quarks (the one corresponding to kq = 0.119). 



To get the final result, that can be confronted to the results of other approaches, we 
need to convert our values to the MS scheme and express it in the physical units by using 



[m K * 



2.7(1) GeV. After recalling that m* gi • c MS (/z) 



m 



MS 



(//), where the function 



c (/•*) is known to 4-loop accuracy El (see App. of the present letter), one can easily 



solve that equation to obtain the standard value m 



MS, 



m r 



Our results are: 



Lattice — I 



m 



MS, 



m 



MS, 



m c )vwi = 1-144(3) GeV , 
m c ) AWI = 1.373(34) GeV , 



(14) 



Lattice — II 



m 



MS, 



m 



MS, 



m c )vwi = 1-132(3) GeV , 
m c ) AWI = 1.325(43) GeV . 



(15) 



The above results are obtained by using tif = and Aq^^ = 0.25 GeV. We see that for 
both sets of our lattice data the two equivalent methods (VWI and AWI) give different 
results. The reason for that discrepancy most probably comes from the lattice artifacts 
which are oc (arriQ) n (n > 2). One way of seeing that is to include the higher order effects 

25] (see also the discussion in p6fl ). 



at tree-level by employing the so-called EKLM factors 

The leading effect of the EKLM factors to our result is oc (arriQ) 2 and it modifies the RGI 

charm quark mass as follows 



m r c gi [VWI] 



m r c gi [AWI] 



m, 



(VWIK2' 



12 



M VWI) ) 2 



mf [VWI] 



m r c gi [AWI] , 



where m 



(VWI) 



stands for the bare lattice charm quark mass whose value is m 



0.375(25). When the above modification is included and we pass onto m 



MS, 



m c 



(16) 

(VWI) 

from 



our Lattice-I simulation, we get 



mf s (m c )vwi = 1.161(3) GeV, 

mp(m c ) A wi = 1.342(33) GeV . (17) 

After comparing these to the results (H), we see that the inclusion of the tree level 
O ((amo) 2 ) effects makes our two results getting closer to each other, although it is not 
sufficient to remove the bulk of O (a 2 ) corrections. 

2.4 Systematic uncertainties and our final result 

We will now briefly summarize the sources of systematic uncertainties and comment each 
one of them. 

• The most important source of the systematic error are the lattice artifacts of 0((am) n ) 
(n > 2). It is therefore highly important to repeat our calculation on the lattice with 
a smaller lattice spacing. As our central result we will quote the average of the two 



methods (VWI and AWI), for both our sets of data as given in eqs. (14,15). The larger 
statistical error will be attributed to our final result while the difference between any 
of the two methods and the averaged one is included in the systematic uncertainty. 
In other words we have: 



Lattice - I mf s (m c ) = 1.26(3)(11) GeV , 

Lattice - II m^(m c ) = 1.23(4)(10) GeV. (18) 

The results of our two simulations are very close to each other and for our final 
estimate of the charm quark mass we will quote the value obtained from the first 
lattice, which has the larger temporal extension. 

In our analysis, we used the value of the lattice spacing a' 1 = 2.71(11) GeV, as 
obtained from the comparison of the physical rnx* and the one we computed on the 
lattice. We checked that if, instead of the above value, we use a" 1 = 2.83(15) as 
obtained from the fx decay constant, the final value of the charm quark gets larger 
by 2%. 

To be conservative, when taking the average of the masses obtained by using the two 
equivalent methods, we quoted the larger statistical error. However, the final result 
does not contain —1.5% of the error on the renormalization constant Z m which we 
discussed in the text (see eq. (0)). Although the value of the renormalization constant 
Z^f does not suffer from the same uncertainty, to be on the safe side, we will add 
— 1.5% of error to our final result. 

Whenever needed, we used Aq^^ = 0.25 GeV. Varying this quantity by 10% allows 
one to cover all the presently available lattice estimates for Aqqd [pO] , pl|] . The impact 
of that variation on the final charm quark mass value is ±1.6%. 

9 



In ref. JT^], the authors also estimate the discretization errors on the improvement 



coefficient 6 a — bp. We have varied this coefficient by ±0.03 which introduces the 
change in the central value for the charm quark mass by less that 1% (more precisely 
by ±0.6%). 

• As mentioned in the previous subsection, the result obtained from the naive (linear) 
interpolation practically coincides with the one that we get when employing the inter- 
polation motivated by the heavy quark expansion. Since we have the third (heavier) 
quark mass, we checked that the quadratic terms to both formulae (i,ii), or the use 
of the "hybrid" formula (hi) affects the value of the charm quark mass by no more 
than 1% (see tab. [J). 

• For determination of the charm quark mass we used the masses of the heavy-light 
pseudoscalar mesons with the light quark mass (linearly) interpolated to the strange 
one. We checked that the quadratic interpolation in the light quark mass to get Mp s , 
does not make any influence on our final charm quark mass value. In addition, we 
verified that our results remain unchanged if we use the vector heavy-light mesons 
instead of the pseudoscalars. 

We sum the above sources of errors quadratically and obtain 

mJP(m c ) = 1.26(3)(12) GeV , (19) 

which is our final result. We note also that in passing from the RGI to the MS value, we 
might have as well used n F = 4 and the value a s {m T ) = 0.334(22) |27]]. This, however, does 
not modify the above result, namely we obtain m^{mc) = 1.27(3)(11) GeV. 

3 Summary and Perspectives 

In this letter we have computed the charm quark mass on the lattice by taking the advantage 
of the full improvement of the Wilson action and operators by which all the lattice artifacts 
linear in lattice spacing are absent. Therefore, the computation of the charm quark mass 
on a reasonably fine grained lattice is expected to lead to the result close to the continuum 
limit. Results of our simulations, obtained at single value of the lattice spacing, indicate 
that the lattice artifacts 0(a 2 ) are sensible. As a consequence the two equivalent methods 
to compute the quark mass on the lattice yield different values. This fact largely dominates 
the systematic uncertainty of our calculation, and for the precise determination of the 
charm quark mass it is therefore important to go to ever smaller lattice spacings. Finner 
lattices, in turn, require larger lattices (to keep the same physical volume) and thus more 
powerful computing resources. By using APE1000 we plan to make such a study in the 
near future. We stress also that our computation has been performed in the quenched 
approximation. The (unknown) uncertainty introduced by quenching is not included in our 
systematic error estimate. A naive guess points towards a decrease of the quenched value 
by ~ 5%. f\ However we prefer to quote our result as the quenched one and wait for the 
partially unquenched computations to assess the amount of this source of systematic errors. 

5 This "guesstimate" is based on the present observations according to which the unquenching of the 
light quark masses reduces their values by ~ 10%, whereas nib gets smaller by less than 2% 0j. 

10 



Appendix 

The RGI quark mass is defined as 



rgi m™(v) m RI/MOM M r „ m 

m C m (fi) C RVMOM^) > <■ U > 

where the functions c MS (/z) and c RI/,MOM (/z) are known to 4- loop accuracy |L8], ^3J. For 
n F = 0, they read 

c m {fi) = a s {n) i/u (l + 0.68733 a s (/i) + 1.51211 a s (» 2 + 4.05787 a,^) 3 ), (21) 

c RI/MOM (/i) =a s (/i) 4/11 M + 2.02067 a s (/x) + 14.21925 a s {^) 2 + 138.30689 a s (^) 3 J (22) 
where, for short, we write a s (fi) = a s (fi)/ir. 
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